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A LOWER BOUND ON THE COSMIC BARYON DENSITY 

David H. Weinberg 1 , Jordi Miralda-Escude 2 , Lars Hernquist 3, , Neal Katz 5 

ABSTRACT 

We derive lower bounds on the cosmic baryon density from the requirement 
that the high-redshift intergalactic medium (IGM) contain enough neutral hydrogen 
to produce the observed Lya absorption in quasar spectra. The key theoretical 
assumption that leads to these analytic bounds is that absorbing structures are 
no more extended in redshift space than they are in real space. This assumption 
might not hold if Lya clouds are highly overdense and thermally broadened, but it is 
likely to hold in the gravitational instability picture for the Lya forest suggested by 
cosmological simulations, independently of the details of the cosmological model. The 
other ingredients that enter these bounds are an estimate of (or lower limit to) the 
intensity of the photoionizing UV background from quasars, a temperature T 10 4 K 
for the “warm” photoionized IGM that produces most of the Lya absorption, a 
value of the Hubble constant, and observational estimates of the mean Lya flux 
decrement D or, for a more restrictive bound, the distribution function P(r) of Lya 
optical depths. With plausible estimates of the quasar UV background and D , the 
mean decrement bound implies a baryon density parameter Qf, ^ 0.0125 h~ 2 , where 
h = Ho/100 km s -1 Mpc -1 . With conservative values of the UV background intensity 
and D , the bound weakens to ^ 0.005 h~ 2 , but the clustering of the absorbing 
gas that is required in order to reproduce the observed mean decrement with this 
baryon fraction is incompatible with other properties of quasar absorption spectra. 

A recent observational determination of P(t) implies ^ 0.0125U 2 even for a 
conservative estimate of the quasar UV background, and f^ 0.018h -2 for a more 
reasonable estimate. These bounds are consistent with recent low estimates of the 
primordial deuterium-to-hydrogen ratio (D/H)p, which imply ~ 0.025 h ~ 2 when 
combined with standard big bang nucleosynthesis. Since the bounds account only 
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for baryons in the warm IGM, our results support earlier claims that this component 
is the dominant reservoir of baryons in the high-redshift universe. The P(t) bound 
on fifc is incompatible with some recent high estimates of (D/H)p unless one drops 
the assumptions of standard big bang nucleosynthesis or abandons the gravitational 
instability picture for the origin of the Lya forest. 


Subject headings: intergalactic medium, quasars: absorption lines, cosmology: theory 


Introduction 


Following the discovery of the first z > 2 quasar ( Schmidt 1965| ), Gunn & Peterson (1965) 
derived a stringent upper bound on the density of uniformly distributed, neutral hydrogen in 
intergalactic space, by showing that the redshifted Lya absorption of neutral gas with more than 
~ 10 -4 of closure density would turn quasar spectra virtually black at short wavelengths, contrary 
to observation. They concluded that the intergalactic medium (IGM) must be highly ionized or 
extremely rarefied. Within a few years, it became clear that the ubiquitous absorption lines in 
quasar spectra are predominantly those of intervening neutral hydrogen (Lynds 1971; Sargent et 
al. 1980| ), and subsequent studies have shown that these lines significantly depress the mean flux 
received from high-redshift quasars blue-ward of the Lya emission line ( pke fc Korcyansky 1982 


Steidel fe Sargent 1987|; |Jenkins &: Ostriker 1991| ; Press, Rybicki, & Schneider 1993, hereafter |PRS| ; 
Zuo Sz Lu 1993 ; Dobrzycki fe Bechtold 1 99?] ; Rauch et al. 1997| ) . Furthermore, it is now recognized 


that the ambient ultraviolet (UV) radiation background produced by high redshift quasars will 
strongly photoionize gas near the cosmic mean density, so that a small amount of diffuse neutral 
hydrogen corresponds to a much larger amount of total hydrogen (e.g., Haardt & Madau 1996, 


hereafter HM). In this paper, we will argue that matching the observed Lya absorption leads to 
interesting lower bounds on the mean baryon density of the universe, which can be derived from 
quite general assumptions about the state of the absorbing gas. 

Recent cosmological simulations suggest that “Lya forest” lines arise in diffuse, but 
non-uniform, intergalactic gas, and that they therefore represent a phenomenon closely akin to 
the “Gunn-Peterson effect” ( Cen et al. 1994 ; Petitjean, Miicket, &: Kates 1995 ; ^hang, Anninos,~fc 
Norman 1995| ; Hernquist et al. 1996| ; for related semi-analytic modeling see Bi 199S ; Bi fe Davidsen 


1997: Hui, Gnedin, &; Zhang 1997). Quantitative analyses show that these simulations require a 


high baryon density in order to reproduce the observed mean opacity of the forest (Hernquist et 


ah 1996 ; Miralda-Escude et al. 1996 ; |Croft et al. 1997 ; Rauch et al. 1997 ; ^hang et al. 199?| ). For 
the UV background predicted by HM based on the observed population of quasars, matching the 
mean opacity estimates of PRS typically requires B ^ 2, where 


B = 


n b h 2 

0.0125 “ 3.4 x 10” 10 


V 


(1) 
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is the baryon density scaled to the fiducial big bang nucleosynthesis estimate of Walker et al. 
(1991). Here h = Hq/ 100 km s _1 Mpc -1 , and rj is the baryon-to-photon ratio. 


The bounds on B derived in this paper will not be as high as those derived from the 
simulations, but they have broader applicability because they are not tied to a specific 
cosmological scenario, and the simplicity of the arguments that lead to them makes it easier to see 
how changes in the theoretical or observational inputs affect the final result. We do appeal to the 
simulations to motivate our one crucial assumption, that structures with a volume filling factor / 
in real space have, on average, a filling factor no larger than / in redshift space. This assumption 
can also be phrased as a requirement that typical Lya forest absorbers satisfy X > 1, where 

X _ real space extent H(z)d 

redshift space extent Av 

is the ratio of the Hubble flow across an absorber (with line-of-sight extent d) to its line width Av. 
This assumption would not hold in a model of spatially compact Lya clouds whose line widths 
are determined by thermal broadening. However, in the physical picture that emerges from the 
simulations, the marginally saturated (IVhi ~ 10 13 — 10 13 cnr 2 ) absorption lines that dominate the 
overall absorption usually arise in moderate overdensity structures that are expanding in proper 
coordinates but contracting in comoving coordinates ( Miralda-Escude et al. 1996| , 1997). The 
absorption line widths are determined largely by these coherent internal motions rather than by 
thermal motions, and the velocity extent of these features is generally no larger than the Hubble 
flow across them. 


Empirical support for the assumption that X > 1 comes from observations of quasar pairs, 
which yield a typical transverse coherence scale ~ 150 h~ l kpc for Lya forest systems at z ~ 2 
( Bechtold et al. 1994 ; Dinshaw et al. 19~94| . 1995). For a typical line width Av ~ 25 km s _1 ( |Hu| 
ct al. 1995| ), this lengthscale would imply X ~ 3/R (for Q = 1, z = 2), where R is the ratio of 
transverse extent to line-of-sight extent in real space. If the absorbers are non-spherical, then they 
are more often intercepted when they are closer to “face-on,” so R might reasonably exceed one 
on average. However, the absorbers would have to be highly flattened, coherent sheets in order 
to reproduce both the observed transverse coherence and the observed line widths while having 
an average X significantly smaller than one. The structures giving rise to the Lya forest in the 
simulations are commonly filamentary, with the transverse coherence scale corresponding to the 
thickness of the filaments. 


Rauch & Haehnelt (1995) used the large transverse coherence scale to argue that the Lya 
forest must contain a substantial fraction of all baryons in the universe at high redshift. Our 
arguments here are different from those of Rauch & Haehnelt - in particular, we work directly 
from observed optical depths instead of from a derived HI column density distribution — but the 
spirit is similar. In §2, we consider the lower bound on B that can be obtained from the mean Lya 
flux decrement ( pke fe Korcyansky 198^ ) alone. In §3 we derive a more restrictive lower bound 
from the distribution of flux decrements (or equivalent optical depths), recently measured from a 
set of seven Keck HIRES spectra by Rauch et al. (1997). Our analytic approach complements the 
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direct comparison to simulations carried out by Rauch et al. (1997), which leads to stronger but 
less general bounds on the baryon density. We discuss implications of our results in §4. 


2. A lower bound from the mean flux decrement 


A uniform IGM with neutral hydrogen density nni produces a Lya optical depth 

„2 


7re“ 


Tu = 


m e c 


faKH (z)n H i, 


( 3 ) 


where f a = 0.416 is the Lya oscillator strength and A a = 1216A is the transition wavelength 
( Gunn fe Peterson 1965|) . The Hubble parameter at redshift z is 


H(z) — Hq Hq(1 + z ) 3 + (1 — Hq ~ Aq)(1 + z ) 2 + Aq 


1/2 


( 4 ) 


where Ao is the cosmological constant A divided by 3Hg. For realistic assumptions about the UV 
background, the IGM is highly photoionized, and the neutral hydrogen density is 


n H n e a(T) 
n Hi =--- 


1.16n^a(T) 

f 


( 5 ) 


where a(T) is the recombination coefficient at the gas temperature T, T is the photoionization 
rate, and nn is the total hydrogen density. Condition (||) enforces balance between destruction of 
HI by photoionization and creation by recombination. In gas with T ^ 10 5 K, collisional ionization 
enhances the destruction rate and lowers nni- The mean value of nn is 


n H = 1-07 x HT 7 (1 + zfB cm -3 , ( 6 ) 

with B as defined in equation (jlj). Equations (|j) and (||) assume a hydrogen mass fraction 
X = 0.76 and a helium mass fraction Y = 0.24. For gas at temperature T 4 = T/(10 4 K) ~ 1, the 
recombination coefficient is 

a(T) = 4.2 x 10" 13 T 4 ~ a 7 cm 3 s _1 , (7) 

. Combining equations (|j) ( ^) yields 

t u = 2.31 x 10~ 4 (1 + ^) 5 (1 + H 0 ^)“ 1 / 2 /i- 1 T 4 -°' 7 r: 4 2 H 2 , ( 8 ) 

where T_i 2 = T/(10 -12 s" 1 ) and we have assumed A = 0 to compute H(z). Equation @ agrees 
with, e.g., equation (36) of HM. 

The mean Lya flux decrement produced by this uniform medium is ZA = (1 — e~ T ) = 

1 — e~ Tu = D u . If the medium is optically thin (i.e., t u <C 1), then clumping the gas tends to 
increase D because the mean neutral fraction at fixed temperature increases in proportion to 
( n H)/( n H) 2 - However, once the gas clumps produce absorption lines with optical depths ^ 1, then 
D is decreased by further clumping, because more neutral atoms are added to saturated regions 


(Abel et al. 1997 
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(where they cannot contribute to increasing D ), at the expense of the interclump medium, where 
the resulting absorption must decrease as matter is moved into the clumps. 

Consider an idealized case in which all the gas collects into non-overlapping, uniform density 
clumps, not necessarily spherical, which have a volume filling factor / and, thus, an overdensity 
1 //. If a clump has the same extent in redshift space as it does in real space, then it produces an 
absorption line with optical depth r c = t u / f 2 and flux decrement D c = 1 — e~ Tc . Since these lines 
fill a fraction / of the spectrum, the mean flux decrement is 

D = fD c = f(l- e~ T ^f 2 ) . (9) 

In the optically thin limit, t u / f 2 <C 1, a Taylor expansion of equation (|9|) yields D = D u /f. In 
the saturated line limit, r u /f 2 3> 1, equation (!) gives D = f. Thus, if t u <C 1, clumping the gas 
(lowering /) increases the absorption initially but decreases it once lines become saturated, as 
expected from our argument above. 

It is straightforward to show (by setting dD/df = 0) that the maximum value of D in 
equation ([9]) occurs when the optical depth of the clumps is T u /f 2 = t c = 1.25643...; f c is the 
solution to the equation t c = \{e Tc — 1). If t u > t c , then clumping can only decrease the overall 
absorption, and D u is the maximum value of the mean flux decrement. If t u < t c , then the 
maximum decrement occurs for filling factor 

/ = (r u /f c ) 1/2 = 0.89214ry 2 , (10) 


implying a maximum decrement 

Ana* = / (l - e- fc ) = 0.71533/ = 0.63818ry 2 . (11) 

This maximum decrement applies to any model with an arbitrary distribution function of the gas 
density, because the general case can be treated as a superposition of cases where all the gas is 
in regions of constant density. Thus, the average absorption cannot exceed that of this optimal 
uniform clump case. 

For an observed value of D, equations (|j) and ([fl]) can be combined to give a lower bound on 
B: 


Brr\ ill - 


65.8 ( 0 § I g) (1 + z)~ 5 / 2 ( 1 + n 0 z) 1 / 4 /i 1 / 2 T 4 °- 35 ri / 1 2 2, for D < 0.71533, ( 12 ) 


65.8 


—ln(l - D) 1 V2 (1 + z )~ 5 / 2 (1 + Qqz) l / 4 h l/2 Tf- :i5 T ] Jy 2 , for D > 0.71533. 


The latter equation applies when D > 1 — e~ Tc , so that the uniform IGM is the optimal case. 
Equation ( |l2[ ) gives a lower bound on the baryon density on the assumption that absorbing regions 
have, on average, the same filling factor in redshift space as in real space. 


If one assumes more generally that the absorbing regions have a constant ratio X of real 
space extent to redshift space extent, then the optical depth of the clumps changes to r c = r u Xj f 2 
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and the spectral filling factor to f /X (or to unity if / > X). One can carry through the same 
reasoning as above to find that the bounds in equation ( |l2| ) simply change by the factor X 1 / 2 . 
As discussed in §1, cosmological simulations and quasar pair observations suggest that X ^ 1 for 
systems that dominate the mean absorption, making equation ([l^) conservative. Models with 
thermally broadened or expanding clouds could have X < 1 and thereby escape the bound (12), 
but the absorbers must be inflated by a factor of four in redshift space in order to weaken the 
bound by a factor of two. 

In order to obtain numerical values of B m j n , we must adopt values of the mean flux decrement 
D , the photoionization rate r„ 12 , the gas temperature T 4 , and the cosmological parameters fio 
and h. PRS, using the high-redshift quasar data of Schneider, Schmidt, & Gunn (1994), find 
D(z ) = 1 — e~ TeB ( z ' ) with r e s(z) = 0.0037(1 + z) 3,46 . The absorption data used to derive this fit 
cover the redshift range 2.5 < z < 4.2. When extrapolated to z = 2, this formula yields D = 0.15, 
in agreement with the value derived by Rauch et al. (1997) from Keck HIRES spectra. Based on 
emission from the observed population of quasars and reprocessing by the Lya forest, HM find 
T_i 2 = 1.4 — 1.6 for 2 < z < 3, and this value agrees well with recent estimates from the Lya 
forest proximity effect ( piallongo et al. 1996| ). In cosmological simulations with a photoionizing 
background, Katz, Weinberg, &; Hernquist (1996) find T 4 ~ 0.6 for gas at the cosmic mean density. 
The solid line in Figure 1 shows the lower bound B m i n using the PRS formula for D(z), r_ 42 = 1-4, 
T 4 = 0.6, h = 0.65, and Ho = 0.3. The derived constraint, B ^ 1, is only weakly dependent on 
redshift in the range 2 < z < 3 because the observed dependence of D on z is canceled by the 


redshift factors in equation (12). Our choice of temperature is conservatively low, since overdense 
gas is usually hotter than gas at the mean density, and since the Katz et al. (1996) simulations 
do not incorporate heat injection during reionization ( Miralda-Escude fe Rees 1994| ). Including 
reionization heating, whose magnitude is theoretically uncertain, could raise the gas temperatures 
enough to increase B m ; n by 20 — 50%. 


Of the various inputs to the lower bound in equation ( 12 ), B m \ n is most sensitive to the 
mean flux decrement D itself. As already mentioned, the value obtained by Rauch et al. (1997) 
at z = 2 agrees almost perfectly with the PRS formula that we have used. At z = 3 the PRS 
formula predicts D = 0.36, and Rauch et al. (1997) measure D = 0.32, which rises to D = 0.35 
after including a theoretically estimated correction for continuum fitting bias. Steidel & Sargent 
(1987) find D = 0.24 for a sample of seven quasars with a mean absorption redshift of 2.64; the 
PRS formula yields D = 0.28 at this redshift. These three studies thus agree to 15% or better. 
However, Zuo & Lu (1993) and Dobrzycki & Bechtold (1996) find values of D that are typically 
about 35% smaller than the PRS values. The Rauch et al. (1997) values seem the most secure 
because they are based on the highest quality quasar spectra, but the sample used is relatively 
small, so further analyses of spectra with similar resolution and signal-to-noise ratio will be needed 
to arrive at a definitive measurement. 


The other significant uncertainty in B min is the value of T_ 42 . The HM estimate is more 
likely to be too low than too high, since it assumes that quasars are the only source of the UV 
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Fig. 1.— Lower bounds on the scaled baryon density B = fl{,/r 2 /0.0125, obtained from the mean 
flux decrement (solid and dotted lines, eq. [jl2[ ) and from the optical depth distribution (filled and 
open circles, eq. |22|). The solid line shows the mean decrement bound with the PRS formula for 
D(z), n 0 = 0.3, h = 0.65, T 4 = 0.6, and T_i 2 = 1.4. The dotted line shows the mean decrement 
bound for conservative parameter choices: D(z ) equal to 65% of the PRS values, h = 0.5 and 
r_i 2 = 0.7. Filled circles show the P(t) bounds at z = 2 and z = 3, using (5 = 0.633, 1% = 0.3, 
h = 0.65, T 4 = 0.6, and T_i 2 = 1.4. Open circles show the corresponding bounds for h = 0.5 and 
r_i 2 = 0.7. Crosses show Rauch et al.’s (1997) estimates of B for h = 0.5 and r_i 2 = 0.7, obtained 
by comparing the observed optical depth distribution to that derived from a numerical simulation 
of the standard cold dark matter model. 




background. However, HM do include an extrapolation of the quasar luminosity function to allow 
for faint sources below existing survey limits, and Rauch et al. (1997) conclude that T _ 4 2 could be 
up to a factor of two below HM’s estimate, at the price of worsening the agreement with proximity 
effect estimates (e.g., Bajtlik, Duncan, fe Ostriker 1988 ; Bcchtold 1994 ; |Giallongo et al. 1991 ; 
Cooke, Espey, Carswell 1997 ). The dotted line in Figure 1 shows the value of B m i n obtained 
after reducing D(z) to 65% of the PRS values, lowering T _ 4 2 to 0.7, and, for good measure, 
dropping h to 0.5. With all of the parameters in equation (Jl2j) pushed to these favorable values, 
the lower bound is B m ; n ~ 0.45. Of course, achieving the observed absorption with B = B m\n 
requires that all of the baryons be collected into uniform density, 6000 degree gas clumps that 
have r = f c , which further requires that the overdensity of these clumps evolve with redshift as 
1// = 0.71533/D(z) (see equation [pd[|). 


3. A lower bound from the optical depth distribution 


The clustering pattern required to achieve the bound of equation (12) is not just physically 
contrived, it is inconsistent with the observed properties of quasar spectra. These show absorption 
at a range of optical depths, not a set of r = 1.25 lines separated by absorption-free regions. The 
two cosmological simulations discussed in Rauch et al. (1997) produce spectra whose optical depth 
distributions agree well with observations, and Rauch et al. argue that they therefore provide 
reliable models of the IGM with which to derive B given values of T_i 2 and h. Here we extend the 
argument of §2 to derive a lower bound on B directly from the observed optical depth distribution, 
based on some simple assumptions about the IGM. 

If we ignore peculiar velocities and thermal broadening, then the optical depth of gas with 
overdensity (uh/uh) and temperature T is 


T = T u 


»H 

n H 


a(T) 

a(T u )' 


(13) 


where t u and T u are the optical depth and temperature that the gas would have if it were 
uniformly distributed. The recombination coefficient a(T) oc T~ 0 -' for the temperature range 
of interest. If we assume a relation between density and temperature, then we can associate a 
density with each optical depth r. In cosmological simulations that adopt B = 1 and the HM 
background spectrum, Weinberg, Hernquist, & Katz (1997) find T 4 « 0.6 {pb/Pb) 0M for the gas 
that dominates the absorption, with higher temperatures in collapsed, shock heated regions. The 
T oc p o e relation arises because denser gas absorbs energy from the photoionizing background 
more rapidly (see further discussion by Croft et al. 1997 ; Hui fc Gnedin 1997 ; Miralda-Escudcj 
et al. 1997). Reionization heating (Miralda-Escude & Rees 1994) can raise the temperature of 
the gas independent of density and thus tends to produce a higher multiplicative constant and a 
weaker trend with density. We will assume 


T 4 = T 4 (n H /n H ) 7 


(14) 
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which implies 

rn = n H ( j/tu) 13 , where (3 = (2 - O. 77 )” 1 . (15) 

Theoretically plausible ranges at z = 2 — 3 are T 4 ~ 0.5 — 1.5 and 7 o.O — 0.6 (^/3 r ^ j 0.5 - 0.65; 


see 


Hui fc Gnedin 1997). Since shocked gas is heated to higher temperatures and produces 


less absorption, assuming that all of the gas lies on the temperature-density relation tends to 
underestimate the amount of baryons present. 


Let P(t) denote the optical depth probability distribution, so that P(r)dr is the probability 
that a randomly selected point on a spectrum has optical depth in the range r —>• r + dr. We can 
use equation (|H]) to compute the mean density by integrating over P(t): 


n H 



n H (r)P(r)dr = n H r u 13 



T ,3 P(T)dT, 


(16) 


and thus 


T u = 


poo 

/ T^P(r)d7 


ii//3 


(17) 


IPo J 

where t u depends on B and other parameters as indicated in equation (pj). We have relied on 
our assumption of Hubble flow broadening both to associate a physical density with a spectral 
optical depth and to identify a probability P(r)dr in redshift space with an equal filling factor in 
real space. We can crudely incorporate a more general model by assuming that the ratio X of 
real space extent to redshift space extent is a function of r alone, X = X(t). If X(t ) > 1, then 
“squeezing” in redshift space has enhanced the optical depth by a factor X(t), but the spectral 
filling factor P(r)dr corresponds to a larger real space filling factor X(r)P(r)dr. Thus, 


poo poo 

nn= n n (r)P(T)dT = n u r~^ [t/X(t )}' 3 X(r)P(r)dr, 
Jo Jo 


(18) 


implying 


T n = 


IJ 0 


) 

[X(r)] l -^T 0 P{T)dT 


It is conceptually helpful to rewrite this equation in the form 


= x (l-(3)/(J 

*■ 1 /1 


IJ 0 


r^P{r)di 


1/0 


1/0 


where X w is the value of X weighted by its contribution to the integral (19): 


X w = 


/ 0 oq [X(t)] 1 -/VP(t)cP 
fo°° T^P(r)dr 


1 i/(i-0) 


(19) 


( 20 ) 


( 21 ) 


In the picture suggested by cosmological simulations, X(r) is typically smaller than unity for small 
r, since the lowest optical depths arise primarily in underdense regions that are expanding faster 
than the Hubble flow. It then rises above unity for intermediate r, where the gas is expanding 
slower than Hubble flow as it falls into overdense structures. Saturated absorption lines are often 
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thermally broadened even in the cosmological simulations, so X(r) may fall below unity again 
at high r. Averaging over these three regimes, X w ~ 1 is a fair approximation for gravitational 
instability models of the Lya forest (see below). Values of X w <C 1 would require dense clouds 
that were thermally broadened or expanding rapidly relative to Hubble flow. 


We can combine equation 
new lower bound on B, 

B min = 65.8 X[}r^ 


3 ) with equation ( 0 ) for the uniform optical depth to obtain a 


r^P(r)dr] ‘ ' (1 + z)- 5 / 2 (l + D^zf^h 1 / 2 ^ 3 V / 1 2 2 . (22) 

Jo 


If all baryons were in the form of intergalactic gas that followed the temperature-density 
relation (14), then equation (22) would be an estimate of the baryon density rather than a lower 
bound. However, some fraction of the baryons should be in stars, shock heated gas, and (possibly) 
baryonic dark matter. Furthermore, since it is difficult to measure the optical depth once e~ T 
is close to zero, we will in practice have to compute the integral in ([ 22 ] ) with a conservative 
assumption, e.g., that all regions with r > r max ~ 3 (transmission less than 0.05) have r = r max . 
Equation (p5) therefore leads to a lower bound on B rather than an estimate. 


Note that if we adopt the optimal uniform density clump model from §2, then 
P(r) = /<^('f'c) + (l — /)d(0), where 5(x) denotes the Dirac-delta function. With / 3=1/2 , the integral 
in equation ( [22| ) becomes f, which, by equation (|ll|), is t^ 2 D/(1 — e~ Tc ) = (D/0. 63818). With 
X w = 1, we recover the bound ( 0 ) from §2, as expected. The maximum value of (1 — e t )/t 1//2 
occurs when r = f c , so for any other optical depth distribution that has the same mean decrement 
D = / 0 °°(1 — e _T )P(r)(fr, the factor / 0 °° t 1 / 2 P(t)cIt must exceed (D/0.63818), and equation ( |22|) 
yields a more restrictive lower bound on B than equation (^). This derivation from the optical 
depth probability distribution is an alternative route to (|l2| ). 


We have checked that equation (22) with X w = 1 gives a lower bound to B in realistic models 
of the Lya forest by applying it to spectra from Croft et al.’s (1997) simulations of three cold 
dark matter (CDM) cosmological models with B = 1. We measure P(t) from the simulations 
at z = 2 and z = 3, in each case using the the value of T_i 2 required to match the PRS mean 
flux decrement. [See figure 11 of |Croft et al. 1997| for plots of P(t) at z = 2.33.] We then apply 
equation (^2|), assuming that all regions with r > 3 have r = 3 in order to account for the limited 
ability of realistic data to estimate optical depths in saturated regions. The derived lower bound 
is about 70% of the models’ true baryon density in the “standard” (SCDM, H = 1, h = 0.5, 
og = 0.7) and open (OCDM, Ho = 0.4, h = 0.65, as = 0.75) models and about 50% of the true 
baryon density in the COBE-normalized H = 1 model (CCDM, Q = 1, h = 0.5, as = 1.2). The 
bounds are insensitive to the value of r max , changing by at most 5-10% if r max is changed to 2 
or to 4. The baryons “missed” by equation (22) are those that have been shock heated to high 
temperatures, condensed into very high density clumps with r>3 (giving rise to Lyman limit 
and damped Lya systems), or converted into stars. These are a larger fraction of the total baryon 
density in the CCDM model because of its higher mass fluctuation amplitude. For the SCDM 
model, we have estimated X w directly by comparing the derived values of B m [ n to values derived 
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from spectra that include no thermal broadening or peculiar velocity distortions (and thus have 
X w = 1 by construction). We find X w = 1.15 at z = 2 and X w = 1.35 at z = 3, so while infall (in 
comoving coordinates) leads to X w > 1 as expected, assuming X w = 1 only reduces -B m in by 4% 
at z = 2 and 9% at z = 3 in this model. 


Rauch et al. (1997) have measured the cumulative distribution of Lya flux decrements from 
Keck HIRES spectra of seven quasars, whose emission redshifts range from 2.5 to 4.5. They find 
that artificial spectra from the SCDM simulation mentioned above reproduce the observed flux 
decrement distribution quite accurately, once T_i 2 is chosen so that the mean decrement matches 
that derived from the data. The artificial (and real) spectra include noise and are measured 
with a locally estimated continuum level. Since we would like to use the true optical depth 
distribution, free of observational artifacts, to obtain R m i n , we measure P(r) directly from the 
noiseless simulated spectra and use it as a surrogate for the observed P(r) in equation (22). As 
before, we set r = 3 in the regions where r > 3. These high optical depth regions cover 9% of the 
spectrum at z = 3 and 2.5% at z = 2. 


The filled circles in Figure 1 show the bound B m i n obtained from equation ( |2^ ) at z = 2 and 
z = 3 assuming X w = 1, (3 = 0.633 (i.e., 7 = 0.6), Qq = 0.3, h = 0.65, T 4 = 0.6, and r _ i2 = 1.4. 
These assumptions correspond to those used in obtaining the solid line, and, as expected, the 
values of B m ; n derived from the optical depth distribution are higher than those derived from the 
mean flux decrement alone. At z = 3 the bound increases by about 25% to B m i n = 1.35, and at 
z = 2 it increases by nearly 70% to R m i n = 1.5. The open circles show the bounds derived for 
h = 0.5 and T_i 2 = 0.7 (with other parameters unchanged): R m i n = 0.85 at z = 3 and B m \ n = 0.95 
at z = 2. These are substantially stronger than the mean decrement bounds represented by the 
dotted line, which also assume h = 0.5 and T_i 2 = 0.7. However, in going from the solid line to 
the dotted line we reduced the mean flux decrements by 35%, while the filled and open circles are 
both based on the Rauch et al. (1997) measurements of the optical depth distribution. If these 
measurements are accurate (the small sample size being their primary limitation), then the mean 
baryon density cannot be much smaller than B ~ 0.9 unless T_i 2 is smaller than 0.7, h is smaller 
than 0.5, or Lya absorbers are substantially larger in redshift space than in real space so that 
X w < 1. The crosses in Figure 1 show the values of B that Rauch et al. (1997) obtain by requiring 
that the SCDM simulation match the observed P{t) with T_i 2 = 0.7 and h = 0.5. These lie about 
a factor 1.5 above the open circles because the direct comparison to the simulation allows these 
estimates to include the contribution of stars and shock heated gas. 


4. Discussion 


The bound on the baryon density from the mean flux decrement, equation (12), and the bound 


from the optical depth distribution, equation (22), are the principal results of this paper. The 
mean decrement bound leads to B ^ 1 if one takes the most convincing estimates of D and of the 
photoionization rate T_i 2 from quasars. With the Rauch et al. (1997) determination of T’(r), the 
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optical depth distribution bound gives B ^ 1 even if one reduces r_i 2 to half of the best estimate 
value. These constraints imply a baryon-to-photon ratio r/ ^ 3.4 x ICR 10 and, assuming the 
standard model for big bang nucleosynthesis, a corresponding primordial deuterium-to-hydrogen 
ratio (D/H)p ^ 6x 10 -5 . This constraint is consistent with the estimate (D/H) p = 2.3±0.4x 10 -5 
of Tytler, Fan, & Buries (1996) from a high-redshift Lyman limit system, but it is inconsistent 
with the estimates of (D/H)p ^ 10 -4 obtained from other high-redshift Lyman limit systems by 
Songaila et al. (1994), Carswell et al. (1994), and Rugers & Hogan (1996ab). Current observational 
estimates of P(r), T_i 2 , and h can only be reconciled with the high (D/H)p estimates by 
abandoning the gravitational instability picture of the Lya forest and returning to a scenario of 
dense, thermally broadened clouds with X w < 1, or, more radically, by abandoning standard big 
bang nucleosynthesis. 


Given an estimated value of T_i 2 , equation ( f±| ) provides an estimate of the number of 
baryons in the warm (T ~ 10 3 — 10° K), diffuse ( p/~p <( 10) intergalactic medium. With the 
Rauch et al. (1997) P(r), HM’s value T_i 2 = 1.4, and plausible choices for other parameters (see 
Fig. 1 caption), one obtains Pigm ~ 1.4 at z ~ 2 — 3 (filled circles in Fig. 1). Within standard 
big bang nucleosynthesis, one must generously increase the estimated observational errors in the 
primordial 4 He abundance even to accommodate B as large as 2 (see, e.g., Hata et al. 1995|) , and 


measurements of (D/H) in the local interstellar medium (see, e.g., McCullough 1992; Linsky et al. 


1995 ) imply B 2.4. This high value of Pigm therefore suggests that the warm IGM contains 


most of the baryons in the universe at these redshifts, as predicted by the cosmological simulations 
and as argued by Rauch & Haehnelt (1995) on different, but physically related, empirical grounds. 


While most of the high-redshift hydrogen resides in the warm, photoionized IGM, most of 
the neutral hydrogen resides in high column density, damped Lya systems. Because we truncate 
P(t) at r max = 3, equation ( | 22 | ) severely undercounts the gas in these systems, which contribute 
Pdla ~ 0.1 — 0.3 h 2 at z ~ 3 ( Wolfe et al. 1995|) . Our lower limit to Pigm far exceeds the baryon 
density of stars in bright galaxies today, P 0.16 h 2 ( |Persic fe Salucci 1992|) . Thus, the baryons 

that were in the warm IGM at z = 2 must either (a) remain in the IGM today, (b) have formed 
brown dwarfs or another form of baryonic dark matter since z = 2 , or (c) have formed stars in 
systems of very low surface brightness that have been missed in standard estimates of the galaxy 
luminosity function. 


There are a number of anticipated observational developments that might strengthen 
(or weaken) one’s confidence in the bounds plotted in Figure 1. The most important will 
be determinations of D and P(r) from larger samples of high resolution, high signal-to-noise 
ratio quasar spectra, since the values of B m \„ depend primarily on these observational inputs. 
Ongoing quasar surveys and, in a few years, the Sloan Digital Sky Survey, will yield improved 
determinations of the quasar luminosity function, which can be combined with the HM formalism 
to yield more definitive estimates of the quasar contribution to the photoionizing background. 
Further analyses of quasar pairs and studies of absorption line shapes in high resolution spectra 
may provide more compelling evidence for extended Lya forest absorbers that are broadened 
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largely by Hubble flow, as predicted by cosmological simulations. Such observations would support 
our key theoretical assumption that X w ^ 1, and they would reinforce the view that the Lya 
forest arises in a smoothly fluctuating intergalactic medium that is the dominant reservoir of high 
redshift baryons. 


We thank Michael Rauch and Rupert Croft for helpful discussions. We acknowledge support 
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